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Abstract. The Bel and Bel-Robinson tensors were introduced nearly fifty years 
ago in an attempt to generalize to gravitation the energy-momentum tensor of 
electromagnetism. This generalization was successful from the mathematical 
point of view because these tensors share mathematical properties which are 
remarkably similar to those of the energy-momentum tensor of electromagnetism. 
However, the physical role of these tensors in General Relativity has remained 
obscure and no interpretation has achieved wide acceptance. In principle, they 
cannot represent energy and the term superenergy has been coined for the 
hypothetical physical magnitude lying behind them. In this work we try to 
shed light on the true physical meaning of superenergy by following the same 
procedure which enables us to give an interpretation of the electromagnetic energy. 
This procedure consists in performing an orthogonal splitting of the Bel and Bel- 
Robinson tensors and analysing the different parts resulting from the splitting. 
In the electromagnetic case such splitting gives rise to the electromagnetic energy 
density, the Poynting vector and the electromagnetic stress tensor, each of them 
having a precise physical interpretation which is deduced from the dynamical laws 
of electromagnetism (Poynting theorem). The full orthogonal splitting of the Bel 
and Bel-Robinson tensors is more complex but, as expected, similarities with 
electromagnetism are present. Also the covariant divergence of the Bel tensor is 
analogous to the covariant divergence of the electromagnetic energy-momentum 
tensor and the orthogonal splitting of the former is found. The ensuing equations 
are to the superenergy what the Poynting theorem is to electromagnetism. Some 
consequences of these dynamical laws of superenergy are explored, among them 
the possibility of defining superenergy radiative states for the gravitational field. 



PACS numbers: 04.20.Cv, 04.20.-q, 04.40.-b 

1. Introduction 

General Relativity is, in some aspects, a peculiar theory. In it the spacetime itself 
is part of the degrees of freedom and this fact brings to General Relativity some 
complications not present in other theories where the fields are set in a fixed spacetime 
background. One of these complications is the impossibility of defining a local 
invariant concept of gravitational energy density. The accepted argument to sustain 
this assertion relies on the equivalence principle. The consequence of this is that any 
geometric object representing gravitational "energy-momentum" can always be set to 
zero in a suitable coordinate system or frame and this property cannot be fulfilled 
by a tensor. Only a pseudo-tensor can accomplish this task but gravitational energy- 
momentum pseudo-tensors are not unequivocally defined because, by the very nature 
of a pseudo-tensor, they are always tied to a given frame or coordinate system. The use 
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of a pseudo-tensor makes it very difficult to address problems sucli as tlie calculation 
of the gravitational energy radiated by a source. 

Different approaches to the "gravitational energy problem" in General Relativity 
have been provided along the years and no general formalism has emerged (although 
formalisms tailored for particular important cases do exist). One of these approaches 
seeks to enhance the formal similarities between electromagnetism and gravitation in 
order to find a replacement for the missing "gravitational energy-momentum tensor" . 
The idea is to take the electromagnetic energy-momentum tensor and translate it into a 
gravitational counterpart by somehow replacing the Faraday tensor with the Riemann 
tensor in the expression giving the energy-momentum tensor for electromagnetism. 
This translation is by no means straightforward due to the different nature of the 
Riemann and Faraday tensors but it can certainly be accomplished. The result of this 
translation is a four index tensor quadratic in the Riemann tensor which was first found 
by Bel [3] . The Bel tensor has mathematical properties which are remarkably similar 
to the electromagnetic energy-momentum tensor (see theorem 14. II for a summary). An 
important particular case arises if we replace in the definition of the Bel tensor the 
Riemann with the Weyl tensor to give the Bel- Robinson tensor [2j . 

From the above considerations it is clear that the Bel tensor will represent a 
quantity which is different from energy. This new quantity was called "superenergy" by 
Bel and its status in General Relativity has been subject to much debate and no widely 
accepted conclusions have been reached. A simple dimensional analysis shows that in 
geometrized units the physical dimension of superenergy is L^'^ where L represents 
length. Another important property is the tensorial character of superenergy. This 
means that if we work with gravitational superenergy instead of gravitational energy we 
can avoid all the technical complications arising when one works with pseudotensors. 
One of the main goals of this paper is to show what the consequences are of considering 
superenergy as a measurable physical quantity on its own. This means that we are 
not concerned in this work with the possible relationship between superenergy and 
other quantities with dimensions of energy. 

In order to carry out our program we need to find the orthogonal splitting with 
respect to an observer of the Bel tensor (so we will be able to explain what the observer 
obtains when measuring superenergy) and we need to find the variation of the different 
parts of the orthogonal splitting along the observer's path. The outcome of this last 
part is a set of equations which we call the dynamical laws of superenergy and they 
are the most important result of this paper. 

We may examine at this point what the above procedure yields in the case of 
electromagnetism. In this case we are working with a quantity with dimensions of 
energy instead of dimensions of superenergy but this is now of no relevance. The 
different parts resulting from the orthogonal splitting of the electromagnetic energy- 
momentum tensor are the electromagnetic energy density, the Poynting vector and 
the electromagnetic stress-tensor. The utility of each of these parts is explained 
in basic electrodynamics textbooks. The dynamical laws of electromagnetic energy 
are contained in the Poynting theorem and it is through this theorem that the 
electromagnetic energy density and the Poynting vector gain their full physical 
meaning as measurable quantities. The Poynting theorem is nothing less than 
the orthogonal splitting of the covariant divergence of the electromagnetic energy- 
momentum tensor. The parts of this splitting are the variation of the electromagnetic 
energy density and the Poynting vector along the observer's path. The Poynting 
theorem enables us to draw conclusions as important as the characterization of 
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radiative electromagnetic fields or the expression for tlie total force acting on an 
electromagnetic system. 

In General Relativity we may consider the expression for the covariant divergence 
of the Bel tensor as the gravitational counterpart of the covariant divergence of 
the energy momentum tensor of electromagnetism. Therefore if we perform the 
orthogonal splitting of the former we will obtain a set of equations which can be 
regarded as the counterpart of the Poynting theorem. As mentioned before these 
equations are the dynamical laws of superenergy and they are far more complex than 
electromagnetism's Poynting theorem. However, we can still follow the same procedure 
as in electromagnetism to draw some conclusions and, for example, we can decide 
in a covariant way when a gravitational system is radiating superenergy (intrinsic 
superenergy radiative state). This was already attempted by Bel in the late fifties but 
since the full set of dynamical laws of superenergy was not available, Bel's result does 
not apply to cases that are sufficiently general. 

The paper is organized as follows: in section [5] we review the notation and 
the essential concepts of orthogonal splittings. In section [3] we find the orthogonal 
splitting of the covariant divergence of the electromagnetic energy-momentum tensor 
in a general spacetime (theorem 13. 2|) . This is the complete version of the classical 
Poynting theorem and some of its consequences are discussed. In section |4] we present 
the Bel and Bel-Robinson tensors and their essential mathematical properties are 
summarized in theorem 14.11 Section [5] contains the orthogonal splitting of the Bel- 
Robinson tensor and we study the basic mathematical properties of the different parts 
of the orthogonal splitting. Since these parts are expressed in terms of the electric and 
magnetic parts of the Weyl tensor, we can obtain particular canonical forms valid for 
some Petrov types fsubsection l5.ip . Section [5] is devoted to the orthogonal splitting of 
the Bel tensor and section [7] contains the main result of this paper which is theorem 
17.11 This theorem spells out the different parts of the orthogonal decomposition of 
the covariant divergence of the Bel tensor (see equation (|4.17|) ) which as explained 
above are the dynamical laws of superenergy. In section [8] we study the radiation of 
superenergy from a general point of view. To that end the definition of an intrinsic 
superenergy radiative state is put forward (definition 18. 3p . 

The main results of this paper rely on heavy tensor calculations which can only 
be carried out with the aid of a computer algebra system. All the calculations of this 
paper have been undertaken with the computer program xAct [33j. xAct is a suite 
of MATHEMATICA packages devised to perform calculations in General Relativity 
and Differential Geometry. Among the many features of the xAct system we stress its 
ability to canonicalize tensor expressions by means of powerful algorithms based on 
permutation group theory (package xPerm), the excellent implementation of tensor 
calculus (package xTensor) and the possibility of working with frames and tensor 
components (package xCoba). In appendix A we provide further details about how 
xAct has been used in this paper. Currently, no other computer algebra system, either 
free or commercial, has the capabilities to perform the calculations needed in this 
paper. 

2. The orthogonal splitting 

We start by introducing the basic notation and conventions which will be adopted in 
this paper. We shall work in a four dimensional smooth Lorentzian manifold V which 
we will call spacetime. The abstract index notation is followed throughout to denote 
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tensors on V with Latin lowercase letters reserved for the abstract indices. We use 
bold typeface for component indices. Round (square) brackets enclosing indices denote 
index symmetrization (antisymmetrization). Unless otherwise stated all tensors are 
assumed smooth and defined globally on V . The metric tensor is gat and our signature 
convention is (— , +, +, +). This metric is used to raise and lower indices in the usual 
way. Associated with the metric is the volume element which we denote by rjatcd- The 
Levi-Civita connection compatible with gab is the only affine connection Va satisfying 
Vagfcc = and our convention for the curvature tensor of this connection is fixed by 
the Ricci identity 

The Ricci tensor and the scalar curvature are Rm = R^ad ^^^'^ = respectively. 
From these, the Einstein tensor is defined by the familiar formula Gab = Rat — Rgab/2. 
The Lie derivative with respect to any vector field X° is the differential operator £x ■ 
Geometrized units with SttG = c = 1 are used unless otherwise stated. The end of a 
proof is marked with D. 

Specially important for us are unit timelike vector fields. For any such vector 
field, the family of its integral curves defines a timelike congruence or observer set. 
This unit timelike vector field enables us to perform an orthogonal splitting (also called 
3+1 decomposition) of any tensor on V. The orthogonal splitting lies at the basis of 
many studies and formalisms in General Relativity and has been extensively studied 
in the literature but since it will be used in this work many times we now review its 
essentials (good accounts can be found in [171 E2])- Let n° be any vector field with 
ria/i" = — 1 and define the spatial metric hab by 

hab = gab + naHb, habh''^ = hac, h°'^=2,. (2.1) 

The tensor hab has the properties of an orthogonal projector. We shall call a covariant 
tensor Ta^ . a^ spatial with respect to hab if it is invariant under h\ i.e. if 

bi ' ' b„^ c^l-'-'i-m ~ ^ bi---b„i J 

with the obvious generalization for any mixed tensor. This property implies that the 
inner contraction of with Ta^ . a^ (taken on any index) vanishes. Wc introduce 
next the orthogonal projection operator defined by 

Pn{La,...aJ^hX---hXLs,...s^, (2.2) 

where Laj^...a^ is an arbitrary tensor. Clearly -P/i(iai...a„) is a spatial tensor. Another 
definition which we need is the generalized inner contraction of the tensor £ai...am with 
the unit normal which is given by 

n''(iai...a„J = n'*! • • • n'*''L...s,...s2...s#j,_, . 

Here J is an ordered subset of the set of abstract indices {ai . . . am} and the dummies 
{si . . . s#j} are placed in those slots of L indicated by J. Therefore n'^ {Lai...am) 
m — # J free indices given by the complement of J with respect to {ai . . . am}- Using 
the orthogonal projection operator and the generalized inner contraction we find that 
any tensor La-^...a,„ can be written in the following way 

iai...a„ = J2 {~^)*-'njPh{n-'{La,...aJ), (2.3) 

Je^({ai...a,^}) 
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where ^({ai . . . a„i}) is the power set of {ai . . . and 

nj = na^- ■ - Ua^, J = {aq, ...,ap}e ^{{ai . . . a^}). 

The right hand side of p.3|) is caUed the orthogonal splitting of La-^^,,,am with respect 
to the unit normal n'^ (we wiU just speak of orthogonal splitting of a tensor if the 
unit normal is understood) . The orthogonal splitting given by (j2.3p is unique and the 
set of spatial tensors {Phi'i^'^ {L ai . . .am))} contains all the information about L^^ 
Equation (|2.3p is just the traditional calculation of the orthogonal splitting of a tensor 
written in a short form. It is possible to study the orthogonal splitting of a general 
tensor in an alternative way if we regard it as a r-fold form (see |38[ 139] for a precise 
explanation of this). 

A trivial example of orthogonal splitting is that of the metric tensor itself which is 
obtained from the first expression in (|2.ip . Another important example of orthogonal 
splitting which is easily deduced from ()2.3p is 

Vabcd = —nabbed + nbEacd — ncEabd + ndSabc, 

where Sabc is the spatial volume element and is defined by 

^abc — ^ Vdabc- 



2.1. Kinematical quantities 

As we explained above, the set of integral curves of n" represents a family of observers. 
In physical applications it is important to introduce quantities describing the relative 
motion of each curve of the family and this is the role of the kinematical quantities. 
To define them we write down the orthogonal splitting of V afT-b which is 

Va^b = -AbUa + ^0hab + <Jab + ^^ab- (2.4) 

The tensor At is the acceleration, the scalar 9 is the expansion and aab, ^ab are the 
shear and the rotation respectively. From the previous equation it is easy to obtain 
expressions for the kinematical quantities in terms of n° 

Q 

A'' = n"Van^ e = Van", UJab - /ifcfVdnc, (Tab = \tK)"^dnc - -hab- (2.5) 

Straightforward properties of the kinematical quantities are 

"■(ab) = CTab, W[ab] = UJab, Cr'^a = 0' (2-6) 

Sometimes the rotation is replaced by the vorticity which is defined as follows 

UJa = -^^abcOj'^'^ UJab — ^abcOJ^- (2.7) 

Each of the kinematical quantities has a precise interpretation which deals with the 
relative motion of the observers of the congruence (see e.g. [TTj [H] for a more detailed 
description of these concepts). 
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2.2. Cattaneo operator 

Another very important object, which is needed when working with orthogonal 
spUttings is the Cattaneo operator also known as spatial connection [14J. If £ai...a„ is 
a covariant tensor then we define the linear operator 

DaLa,...a^=PhiVaLa,...aJ, (2.8) 

with obvious definitions for contravariant and mixed tensors. The Cattaneo operator 
is not a linear connection on the spacetime manifold V because it does not satisfy 
the Leibnitz rule unless both factors of the product upon which Da acts are spatial. 
From its definition, it is clear that DaTi,-^ i,^ is a spatial tensor. Important additional 
properties of the Cattaneo operator are 

Dahbc = 0, DaDb^ - DbDaip - 2n''ujab^c(p = 2uJab£n(p, ^ € C\V). (2.9) 

The Cattaneo operator enables us to write in a compact form the orthogonal 
splitting of any expression involving derivatives. It is specially important in this work 
to find the orthogonal splitting of the covariant derivative of a spatial tensor. To 
illustrate how this works, let us study a particular simple example. Consider Va^b, 
where Lf, is an arbitrary spatial covector {n°-La = 0). In this case formula (|2.3p yields 

S/aLb = DaLb - naPhinPVpLb) - nbPh{n^^ aLp) + UaUbvPn'^V pLq. (2.10) 

Next we use in this equation the relations 

n^VcLa = £nLa - , n^V aLp = -LpVaU'', (2.11) 

and replace in (|2.10p the covariant derivatives of the unit normal by the expression 
given in (12. 4p . After some manipulations equation (|2.10p becomes 

Vaib = -A^LcUbUa + {^^bO ^ £ nLb + {abc + t^&c)^ ^^a + DaLb + 

+nb (^La9 + L^CTac + C^acU , (2.12) 



which has the form of (|2.3p and hence is the complete orthogonal splitting of VaLb- 
Note that £nLa is a spatial covector if La is spatial, due to the property = 0. 

The procedure followed to obtain ()2.12p is easily generalized for the covariant derivative 
of any spatial tensor (see appendix A for more examples). This kind of calculation is 
extensively used in section [T] 



3. Electromagnetism as a working example 

As a preparation for the study which we are going to undertake of the gravitational 
field, we analyse first the case of electromagnetism. The electromagnetic field is 
described by an antisymmetric rank-2 tensor Fab (the Faraday or electromagnetic 
tensor) which satisfies the Maxwell equations 

V[aFbc]=0, VaF%=Jb, (3.1) 

where is the charge current four vector (here we follow the Heaviside-Lorentz units 
system). A very important object in electromagnetic theory is the energy-momentum 
tensor of the electromagnetic field, given by 

Ta' = \{FadF'"' + F:aF*'"') = FadF"'' - ^S^F^aF'". (3.2) 
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In this formula F*^ is the Hodge dual of Fab defined by 

^afc = -jVabcdF"'^- 



Theorem 3.1 The tensor Tab has the following properties 

(i) Tj^ab) = Tab ■ 

(ii) Tab always satisfies the dominant energy condition, namely, for any pair u"" , 
of causal future- directed vector fields the inequality TabU°'v^ > holds. 

(Hi) If the Maxwell equations hold then we have 

VbT\ = Fa'jb, Vaf = 0. (3.3) 

Roughly speaking, the first equation of p.3p tells us that the variation of the 
electromagnetic energy-momentum equals the work performed by the charge current 
and the second equation is the equation of charge conservation. An adequate 
understanding of these informal assertions can be achieved by finding the orthogonal 
splitting of (|3.3p . As an aside remark, we note that the first equation of (|3.3p is not 
in general equivalent to Maxwell equations as is sometimes wrongly stated. 

To find the orthogonal splitting of (|3.3p we first need to find the orthogonal 
splitting of Fab- Define the spatial tensors 

Ea = Fabn\ Ba = F:,n'' (3.4) 

These are the electric and magnetic parts of the Faraday tensor and they characterize 
it completely. The orthogonal decomposition of the Faraday tensor in terms of Ea 
and Ba reads 

Fab = EbUa - EaUb - BpSab- (3.5) 

Using this expression, we can find the orthogonal splitting of the energy-momentum 
tensor Tab which results in 

Tab = UUaUb -f 2P(anb) + Tab, 

U = ^{EaE'' + BaB''), Pa = SabcB''E', Tab = Uhab - EaEb - BaBb- (3.6) 

Also, the orthogonal splitting of is easily found yielding 

f = pn'' -r, p= -fna, r = -h'^'jc 

where p is the charge density and J° is the spatial charge current. Next we replace 
the decomposition of Tab and j°- in (|3.3p and calculate the orthogonal splitting of the 
resulting equations. To achieve this we need to find the orthogonal splitting of VaU, 
\7aPb and Vq7{,c which is done by using the appropriate generalizations of (j2.12p (see 
the proof of theorem 17. li theorem 17.21 in appendix A and especially equation (jA.ip ). 
The final result is presented next. 

Theorem 3.2 The following set of equations 



VbT\ = Fa" 3b, Vaf = 0, 
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is equivalent to 



£nPa = -SabcB^r + EaR + 2eabcP^0j' - PaO - A\Uhab + Tab) - DbT^a 
£nP=-A^Ja + 0p + Dar. 



(3.7) 



(3.8) 
(3.9) 



Equations ()3.7|) - (I3.8|) are presented in basic electrodynamics books under the heading 
of the Poynting theorem and they reflect the transfer of energy-momentum in a system 
composed of charged particles and electromagnetic fields. Indeed, equations (|3.7p - (|3.8p 
provide the well-known physical interpretation of each of the quantities appearing in 
equation (|3.6p : U is the electromagnetic energy density, P" is the Poynting vector 
and Tab is the stress tensor of the electromagnetic field (see e.g. [26] for detailed 
explanations about the role of each of these quantities). 

We must note at this point that equations (|3.7p - (|3.8p are usually presented under 
the assumption that the spacetime is flat and n"' is chosen in such a way that all 
the kinematical quantities vanish. The resulting equations can be always obtained 
locally in a general spacetime if we recall that we can always construct a vector 
field with the property that all its kinematical quantities vanish at a prescribed 
point (equivalence principle). Therefore, we deduce from these considerations that we 
can classify the terms which appear in p.7p - (|3.8p into two categories: those which 
contain kinematical quantities and those which do not. Terms which do not contain 
kinematical quantities can be regarded as representing intrinsic variations of energy 
or momentum (non-inertial terms) whereas terms affected by kinematical quantities 
can be thought of as depending on the observer n° and we shall call them inertial 
terms in analogy to the inertial forces introduced in the study of accelerated systems in 
Newtonian physics. These considerations, although elementary, will play an important 
role in section [H 

3.1. Coupling of the vorticity and the Poynting vector 

If the vector field n'^ is hypersurface orthogonal then p.7p - (|3.8p assume simpler forms 
which can be found in different places in the literature |42j. The general form of 
(|3.7p is written down in [32] but to the best of our knowledge equation (j3.8p does not 
seem to be present in accessible references. Also some of the consequences of (|3.8p do 
not appear to be widely known. To illustrate this fact, consider the inertial terms in 
p.8p . In ordinary units we find that the left hand side of (|3.8p is the time variation of 
the momentum density and therefore the terms on the right hand side of (|3.8p which 
are coupled to the kinematical quantities can be regarded as inertial forces. Indeed 
equation (|3.8p can be interpreted as an equilibrium condition for an electromagnetic 
system which states that the sum of all (inertial and non inertial) forces acting on the 
system equals zero. 

One of the inertial forces is given by 2eabcP'^^'^ or in three- vector notation 2P x Q 
with "x" representing the vector product. If we consider a gyroscope then we find 
that the vorticity lo'^ is related to the angular velocity of the gyroscope. Therefore we 
deduce that an inertial force exits on an uncharged gyroscope when it is placed in a 
radiative electromagnetic field. In SI units this inertial force can be estimated by 



Eg « —P X (J, 



(3.10) 
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where V is the volume of the gyroscope. Thus we conclude that the flux of 
electromagnetic radiation produces an effect on a gyroscope. We must stress at this 
point that this is an observer dependent effect (as it should be because we are dealing 
with an inertial force) which manifests itself in the fact that the angular velocity of 
the gyroscope depends on the observer. An effect similar to this was pointed out in a 
particular case in and this was latter confirmed in [22] . In the former reference it 
was shown that gyroscopes placed in the spacetime generated by a nonrotating charged 
magnetic dipole would precess. As an explanation of this result it was suggested that 
the Poynting vector could cause a measurable effect on a gyroscope's precession and 
it is conceivable that (|3.10p is related to this effect in some way. 



4. Gravitational equations and the Bel tensor 

We start this section by reviewing the well-known formal analogy which exists between 
electromagnetism and gravitation. In this framework the Riemann tensor Rabcd is 
taken to be as the gravitational counterpart of the Faraday tensor Fab and the role of 
the two Maxwell equations is played by the relations 

^[aRbc]df = 0, V dRbpc — 5p&c, Zefa = ^ eRaf — ^ f Rae (4.1) 

The tensor 3a6c is known as the matter current and can be regarded as the counterpart 
of the charge current four vector j". There is an important difference between 
electromagnetism and gravitation in that in the latter we have an extra set of 
conditions: the Einstein field equations 

Gab = 'lab- (4.2) 

Here the tensor Tat is the energy-momentum tensor of the system and must be 
prescribed independently. Clearly any solution of the Einstein equations will be a 
solution of (|4.ip but the converse need not be true. From (|4.ip we derive the well- 
known relation (see e.g. [TB]) 

VaV°i?dc6p = 'VbZdcp - ^pZdcb - ^R^^^R^^j^^ - Rfj'^Rdcpa + Rdcba^p"^ " "^^cape^db + 
+ '^RbecaRdp'': (4-3) 

which can be shown to be a hyperbolic equation for the Riemann tensor. A result 
due to Lichnerowicz [31] proves that if the Cauchy data of (|4.3I) satisfy (|4.2p then so 
does the solution of the hyperbolic equation. Hence, with the provision imposed by 
the Lichnerowicz result, we can regard (|4.ip and (|4.2p as equivalent. 



^.1. Orthogonal splitting oj the Riemann tensor 

The orthogonal splitting of the Riemann tensor was first studied in [1] and since then 
it has been used in many places. Define the left, right and double dual of Riemann 
tensor in the standard fashion 

J^abcd — ^ lo-bpqJ^ cd^ ^abcd — 2 'P?'^'' a& ' ^abcd ~ 2 'ab ^pqcd- 

Next we introduce the following spatial tensors [4] 

Yac = RabcdTl'ri''^, Zac = * Rabcdfl'' n"^ , Xac = * R*abcd'^^ (4.4) 

The symmetries of Riemann tensor entail the properties 

X^ab) = Xab, Yi^ab) = Yab, Z\ = 0. (4.5) 
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These tensors contain all the information in the Riemann tensor as is easily checked 
by a simple count of their total number of independent components. They also enable 
us to find the orthogonal splitting of the Riemann tensor which reads 

Kbcd = 27icn[aYh]<j + 2ft,^[^X^]j, + 2ndn[iY^^^ + 2n[aZ''^]e^f,^ + 2n[bZ^^]£cde + 

+htd {KcX\ - X^J + V {X,d - KdX\) ■ (4.6) 

From this expression is easy to get the orthogonal splitting of the Ricci tensor which 
is 

Rac - Z'^'^EcdbUa + n^y'dUa - Xac " Yac + U^Z'^^'eadb + hacX^'^. (4.7) 

The Weyl tensor Cabcd has the same algebraic properties as the Riemann tensor and 
in addition it is completely traceless. Therefore to find its orthogonal splitting we 
proceed along the same lines as with the Riemann tensor but using different names 
for the tensors introduced in (|4.4p . The precise correspondences are (in the next 
equation Xab, Yab, Zab are defined as in (|4.4p with the Riemann replaced by the Weyl 
tensor) 

Bab = Zab = Z(^ab)i ^ab = Yab = —Xab, E°'a = 0- (4-8) 

The tensors Eab and Bab are known as the electric and magnetic parts of the Weyl 
tensor and they completely characterize the former. Equation (|4.6p becomes for the 
Weyl tensor 

Cabcd = '^n^r^aE^^ ~ 2/i^[^£;^]j, + 2ndn[i,E^]^ + 2nydB\fabe + "^-n^bB^^Ecde + 2h^^^E^^^. 

(4.9) 

^.2. Orthogonal splitting of the matter current 

The orthogonal splitting of Zabc can be calculated if we insert in the last expression 
of (|4.ip the orthogonal decomposition of the Ricci tensor (|4.7p . In this calculation the 
orthogonal splittings of VaXbc, VaYbc, ^ aZbc, ^ a^bcd must be used (see appendix A 
for the explicit expressions). The result is 

5e/a = -Lfnena+ LeUfUa + Jf^Ua+UfJ^a " neJfa+jefa, (4-10) 

where 

Je/ = 2iX^f + Y[/)(ae], + UJe]a) - 2(X"„ + F^JWe/ + 2eab[eDf]Z-\ (4.11) 
J, a ^ 21; "a,, + X'.aae + ha, (- J^'b^ + ^"V,,^ + (-2X,^ + Ya')<J,, + 

+2Y,'ojab - {X\ + 2Y\)u:a, - Ya'u,,,+ sM-A' Z\ - DaZ'^ + D^Z\) + 

+eabc{-A''Z% + D,Z^-) + ^{Xa,e - i{£nYae)), (4-12) 

Le ^ A,{X\ + YD ~ A'^iX.a + + 2^"^[a.] + ^ abcZ"' o + D^Y^, (4.13) 
jefa = 2W/Z[ea] + 2We^[a/] - 6^'' (/l^/ ^[eb] + ^ae%/]) + 4t^a^[e/] + 

'^^bcyZ'' Q/le]a^ + ^e]a) + 'iK[fD,^X\ + 2DyX,^a + Zi^l/^Ja- (4-14) 

From these expressions we deduce the properties J[ab] — Jab, j[ab]c = jabc- 
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4-3. The Bel and Bel-Robinson tensors 

Finding a gravitational equivalent of the electromagnetic energy-momentum tensor 
Tab proves to be a delicate issue. The reason for this lies in the impossibility of a 
local definition of the gravitational energy-momentum density due to the equivalence 
principle. Therefore it is clear from the very beginning that any tensor qualifying as 
the gravitational counterpart of Tab must represent a physical quantity different from 
"energy-momentum" . If we are unwilling to introduce "new quantities" in physics 
then the point of view traditionally adopted consists in resorting to quantities defined 
non locally or using pseudo-tensors (a very good review of the research carried out 
in this direction is [41]). However, if we are ready to deal with a quantity different 
from "energy-momentum" then we find that it is possible to construct a tensor whose 
mathematical properties are similar to the electromagnetic tensor Tab and this is the 
Bel tensor. 

The Bel tensor was first introduced in _3_ in connection with the construction 
of covariant divergences of quantities quadratic in the Riemann tensor. The original 
definition given by Bel can be shortened to the expression 

Ba'f ^ liKpcR* + *Rapc, *R'"""' + *R:pc, *R* + RapcR''"'), (4.15) 

which is formally similar to the first equation in (j3.2p although with more terms due 
to the fact that the Riemann tensor has two blocks of antisymmetric indices. 
If we expand the duals in (|4.15p we get 

Babcd — RaecfRbd ~ '^^dcRaefpRb'^'^^ ~ l^dbaRcefpRd^^ + Rb J Rdf ae + 

+ \9ba9dcRefphR'^'''' ■ (4.16) 

The Bel tensor has a number of remarkable mathematical properties which are 
summarized next. 

Theorem 4.1 The following statements hold true for the Bel tensor 

(i) Babcd — B(^ab){cd) — Bcdab, ^^Qcrf " 0- 

(ii) (Generalized dominant property) If u^, Uj, u^, are arbitrary causal, future 
directed vectors then BabcdUiU2U3u'l > 0. 

(Hi) Babcd = Rabcd = <J=> 3 fl timelike vector w° such that BabcdU°'u^u'^u'^ = 
0. 

(iv) Equation entails 

1 ^aef 

^aB\^d — Zd^Rbcca + Z,°'^Rbeda — 'jZcdZ'^'^ Rbfae, ^a^bc"^ — 0- (4-17) 

The similarity between the mathematical properties of Bated presented in this theorem 
and those of Tab given by theorem 13.11 is apparent. Therefore the Bel tensor fulfills 
the basic mathematical requirements needed for it to be regarded as the gravitational 
counterpart of the energy-momentum tensor in electromagnetism. In vacuum, the Bel 
tensor acquires a simpler expression which is 

Tabcd = Ca^d ^bpcf + CJ'JCbpdf — -^gabScdCpqrsC^'^^'' , (4-18) 

where Rabcd = Cabcd has been used. The tensor Tabcd is known as the Bel-Robinson 
tensor [2] and it can be defined in any spacetime, whether vacuum or not, by means 
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of equation (|4.18p . All the properties of theorem 14.11 except point are also true for 
the Bel-Robinson tensor with the following changes: Tabcd is totally symmetric and 
trace-free and in point plz|l the Riemann tensor must be replaced by the Weyl tensor. 
If Rab = Agab (Einstein space) then the covariant divergence of the Bel-Robinson 
tensor takes a particularly simple form. 



A full account of the properties reviewed here of the Bel and Bel-Robinson tensors 
together with their proofs can be found in [38] and |10) . In the former reference 
a generalization of (|3.2p and (|4.15p valid for any tensor is put forward. Tensors 
resulting from this generalization are called superenergy tensors and they all fulfill 
the generalized dominant property (generalized dominant superenergy condition). 

What about the physical role of the Bel tensor? This question has been addressed 
many times in the past and no definitive answer exists. Bel himself proposed the name 
of superenergy for the physical quantity which might lie behind the Bel tensor (this 
physical quantity would be represented by the components of Bel tensor in a suitable 
frame). If we denote by L the basic unit in the gcomctrized system then from the 
definition of the Bel tensor we deduce that the physical units of superenergy are 
which can be interpreted as either energy density squared or energy density per unit 
area. Both interpretations have been researched in the literature and the opinion 
favoring the second interpretation seems to have gained weight. For a history of the 
different interpretations of the Bel tensor which have been studied in the past see |38| 
and references therein. 

In the case of electromagnetism we have seen that a full understanding of the 
physical properties of the electromagnetic energy-momentum tensor can be achieved 
by the Poynting theorem. This theorem is nothing but the orthogonal splitting of (|3.3p 
and the different equations of this splitting inform us of the evolution of the different 
parts of the electromagnetic energy-momentum tensor. Therefore it is expected that 
the orthogonal splitting of equation (I4.17|) will yield valuable information about the 
true physical role of the Bel tensor. The calculation of such an orthogonal splitting is 
accomplished in the forthcoming sections. 

5. Orthogonal splitting of the Bel-Robinson tensor 

Before studying the general case of the Bel tensor we calculate the orthogonal splitting 
of the Bel-Robinson tensor. The different parts of the splitting take simpler forms and 
they will give us valuable insights about the general case. To calculate this splitting 
we insert the expression for the orthogonal splitting of the Weyl tensor given by (|4.9p 
into (|4.18p . After some computations we get 

Tabcd = WnaUbUcUd + -iV {anbUcnd) + Gt(^abncnd) + -iQ (abend) + tabcd, (5.1) 



W = EabE'^" + BabB'^", Va = 2Bp'E,ieP\ tab = Whab ' 2{B^-Bbc + E^'Ebc), 



la 



= 0. 



bed 



where 




(5.2) 



t. 



abed — 



^{Bab^cd + T^ab^cd) - hcdtab + '2hb{dtc)a + 2ha(dtc)b " ^ab^crf + 
+ W(habhcd - 2ha(chd)b)- 
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Some of these quantities have been obtained before and have found diverse 
appHcations. The scalar W (superenergy density) and the spatial vector V^, called 
the super-Poynting vector, were first used in [5] to define intrinsic radiation states in 
gravitation theory (see section [5] for more details about this) and the tensor tab was 
used in [11] to show the causal propagation of gravity in vacuum (also the role of tab 
in the definition of radiation states was discussed in this reference). We establish next 
the basic algebraic properties of these quantities 

Proposition 5.1 The following basic algebraic properties hold 

(^) l'{ab) tab: Q (abc) Qabc} l'(abcd) ^abcd: 

(^^) t\ = W> 0, Q\, - Vb, t\,^ = tbc, 

(Hi) Qabc o,nd tabcd Contain all the information about the Bel-Robinson tensor. 

Proof : Points Q and Q can be proven directly from the tensor expressions given 
in (|5.2p but it is far more easier to use (|5.ip and write each part of the decomposition 
in terms of the Bel- Robinson tensor Tabcd- The result is 

W = Tabcdn'^n^n'^n", Ve - -T.bcdh'an'n-n", tab = Tp,„/^P>«,n'■n^ (5.3) 

Qabc — Tpqrsh^a^'^b^ c^ ' tabcd — Tpqj, gh^ ^. (5-4) 

The symmetries expressed in point Q are now a consequence of the total symmetry 
of Tabcd- Point Q is straightforward either from (15. 2|) or from (I5.3|) - (|5.4|l and the 
complete tracelessness of the Bel-Robinson tensor. Thus given tabcd and Qabc if is 
evident from their algebraic properties that we recover the remaining parts of the 
orthogonal decomposition of the Bel- Robinson tensor which proves point (fm| . D 

Remark 5.1 We can obtain an independent proof of point (fm|l of the previous 
proposition if we count the number of independent components of Qabc and tabcd 
and compare their sum with the number of total independent components of Tabcd- 
The respective numbers are 

number of independent components of Tabcd — 25, 
number of independent components of tabcd = 15, 
number of independent components of Qabc — 10, 
10+ 15 = 25. 
Proposition 5.2 

tabcd = ^ iah = ^ T4^ = ^ Cabcd = 0, 

Proof : From proposition 15 . II we deduce tabcd = tab = W — Q. To prove 
the converse, let us assume that W = Q. In this case the first equation of (|5.3|) implies 

Tabcdn'^n^'n^n'' = 0. 

Combining this with point (fmj) of theorem 14.11 applied to the Bel-Robinson tensor we 
deduce Cabcd = 0. Trivially, Cabcd = implies Tabcd = and thus tab, tabcd vanish as 
well. □ 
The importance of this result lies in the fact that evaluation of any of the 
quantities W , tab, tabcd enables a observer represented by the unit timelike vector 
n° to decide if the purely gravitational part of the Riemann tensor (or the Riemann 
tensor itself if we are in a vacuum spacetime) is present or not. Also the variation of 
these quantities along the integral curves of should give a measure of how the Weyl 
tensor changes for this observer. We will turn back to this important point in section 

m 



Dynamical laws of superenergy in General Relativity 



14 



5.1. Canonical forms for the different Petrov types 

We can obtain more interesting properties of the quantities introduced in (j5.2p if 
we set up a suitable orthonormal frame. Such a frame arises in the calculation of the 
canonical forms which Eab and Bab take for the different Petrov types. These canonical 
forms are reviewed in appendix B and we refer the reader to this appendix for more 
details. The results presented in this subsection are algebraic in nature and should be 
understood as formulated in the tangent space of a point. 

Proposition 5.3 The tensor Qabc vanishes if and only if Eab, Bab are linearly 
dependent. 

Proof : From (|5.2p it is easy to show that Qabc is zero if Eab and Bab are linearly 
dependent. Now, if Qabc = then from point (O of proposition 15. II we get 7"* = 0. 
This last condition can be re-written in the form 

Ea'-Brb - E.'-Bra = 0, (5.5) 

from which we conclude that the endomorphisms represented by E'^j^, B\ commute. 
This is only possible for Petrov types I and D as can be easily checked using the 
canonical forms of appendix B (alternatively, two symmetric endomorphisms have 
a common basis of eigenvectors if and only if they commute). For Petrov type D 
trivially Eab and Bab are linearly dependent, so we will assume that the spacetime is 
of Petrov type I. In the orthonormal frame of (jB.ip we find that the only nonvanishing 
component of Qabc is 

Ql2Z = — 2(_Bii_E22 — B22E11), 



and hence (5i23 — implies B11E22 = B22E11 from which we deduce from (jB.ip that 
Eab and Bab are hnearly dependent (recall that -Eii+£'22+-E33 = -B11+-B22+-B33 — 0). 
□ 

From this result we deduce that Qabc resembles in its mathematical properties 
the electromagnetic Poynting vector. We will see later that if we are to study the 
radiation of superenergy then Qabc (or any equivalent tensor thereof) will take over 
the role of the Poynting vector. 

Proposition 15.31 admits the following corollary. 



Corollary 5.1 

(i) Qabc 7^ =^ Petrov type is either II, III, N or I. 
(a) If Petrov type is II, III, or N Qabc 7^ 0. 
(Hi) Petrov type D is the only type in which Qabc always vanishes. 

Proposition 5.4 The following algebraic properties hold 
(i) For Petrov type III we have 

tab — -nabVy — — , Wabc — on-{bc' a) ^ 



. MVaVbVcVd , 12 p^npp--^' 

tabcd — ^ r^n-iabPcPd), Pa T ^, PaP — ■ 

The two independent principal null directions of the Weyl tensor (see e.g. JJW) 
can be calculated explicitly yielding 

= + ^n^W, = + ^n'^W. (5.6) 
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(ii) For Petrov type N we have 

_ VgPb ^ _ VgPbPc _ VgPbPc'Pd ^ , „ -p "P" — T/f/2 

^ab — ' Wabc — ttto ' ''abed — ti^^ : ' a ' a' — ■ 

In this case the only independent principal null direction of the Weyl tensor is 
k" = Wn'' +r''. (5.7) 

Proof : The proof of this result consists in using the canonical forms for Petrov 
types III and N written in appendix B to find canonical forms for tab, V"", Qabc and 
tabcd- These canonical forms lead then to the expressions presented in points ^ and 

We detail next this procedure for each of the Petrov types. 
- Petrov type III: using the frame of (jB.4[) we get 



-Qi33 — -Q122 — Qui — 2(i3^2+-Bi2), Vi — -2{Ef2+Bl2), h2 = ^33 — 2(i?^2+-Bi2). 

with all the other components of Qabc, T-"^ , tab being zero. From these expressions we 
deduce 



Qabc = hbc'Pa + hac'Pb + hab'Pc— „2 ^ t?2 ' ~ ^ (£'^2 +-^12 ) ^afe ~ r,^ f2 " i ^p2 \ ' ^^■'^ 

-^^12+ ^12 ^(,-^^12 + ^12) 

and using point (0) of proposition 15.11 we conclude 



i'l2+-C'12- ^ I l^aP - — ■ 

Replacing this back in (|5.8p we obtain the expressions sought for tab and Qabc- 
Inserting the values just found for tab in the formula for tabcd of (|5.2p yields 

tabcd = ^{BabBcd + EabEcd) H "^'^ -^ha^d^c) + -^{hcd'Pb — h^dPc — hbcVd)- 

Again using the canonical forms of (|B.4p we transform the term A{BabBcd + EabEcd) 
into 

_^^PaPbPcPd ^ _4_^^^^^^^^^^ ^ ^^^^^^ ^ ^^^^^^ ^ hacVb)Vd). 

Combining the last two equations we find the expression for tabcd given in the 
proposition. It is now a simple calculation to check that the vectors and fcf are 
indeed null and that they fulfill the properties 

Tabcdkik^k^ki = 0, T'(jhcdfc2 ^2^2^2 ~ ^ 

which implies that fc" and fcj are the Weyl tensor principal null directions (see [36] p. 
328). 

- Petrov type N: In this case, we obtain in the frame of (jB.sp 
Qui — Vi— —til = — 4(i?22 + -E'22)' 
with the other components vanishing. Hence 

n — 'Pg'Pb'Pc _ VcVd . pt'^ I — 

Wabc — , „ / j-,2 I 7^9 \2 1 ''cd — , , „2 i 7-12 \ ^22 ^22 — 



VaVbrc , _ VcVd ^ r2 ^ p2 _ _ w_ 

16(512 + EI2? ' ^Bl2 + EI2) ^ ^22 + ^22 - 4 - 4 . 
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Similarly, working in the canonical frame we obtain that the only nonvanishing 
component of tabcd is 

iiiii = 4(^22 + ^22)- 

Combining the previous pair of equations the expressions of point ([np follow. Also it 
is a simple matter to check that is null and that Tabcdk°'k^k'^k'^ = 0. D 
An important result of this proposition is that for Petrov types III and N the 
Bel-Robinson tensor is characterized by just two independent quantities which are W 
and V"' and thus we can say that the number of algebraically independent components 
of the Bel-Robinson tensor is two for these Petrov types. This is not true of the other 
Petrov types and therefore some conclusions drawn from considerations involving type 
III and N might not carry over to other Petrov types. An example of this is the 
definition and study of gravitational radiation using the Bel-Robinson tensor where, 
traditionally, a nonvanishing vector for any observer n° has been regarded as an 
intrinsic state of gravitational radiation ^ (see definition lS.ip . We will see in section[7] 
that this condition is not general enough and indeed in certain Petrov type I spacetimes 
we can still speak of an intrinsic state of gravitational radiation with V"' being zero. 

6. Orthogonal splitting of the Bel tensor 

The orthogonal splitting of the Bel tensor is obtained by replacing the expression for 
the Riemann tensor given by (|46)) in (|4T6| with the result 

Babcd = WnanbUcUd + {aribncUd-^ -|- 2n(^aQb)cd + '^'n[dQc)ab + 

-HabncUd + tcdUaUb + An (^ail) {c""- d) + tabcd- (6.1) 

Each of the spatial parts of the Bel tensor is defined as follows 

W = \{XabX''' + YabY'^'') + ZabZ''\ V, = EabciYa'Z'" - X/Z''^), 

ted = hcdW — X^"'Xda — Yc ""Yda — ZacZ'^d ^ Z^Zda, 

Kd = 2A:(/yb)a - XbdY";^ - YbdX\ + hbdi-X^-^Y^^ + Z^-^Z^^ + X^Y"^) - 

^Z\Z^^ — ^h^^jja 
Qbcd = ^cdPh + 2Z°^ (-^cf ^bae + ^c)6a^''e + <^c)aeXb^ + 

-|-2Z"^ {hb(^c{^d)aeX^ ^ — S^^^fXj) ~ X^(^d£c)ba + hcdXj Ebef ~ Xb(^d'^c)ae^ ■ 

The expression for tabcd is a bit long and is omitted (its explicit form is not needed in 
this paper). 

Proposition 6.1 The tensors tab, t*^^^, Qabc '^'^'^ iabcd satisfy the following basic 
algebraic properties 

't(ab) — tab, t(ab) — ^ab, Q a{bc) — Qabc, t{ab)cd = tabcd = tcdab, t ^ — W , Qb a — T^b, 
^ abc — ^bc 

Proof : These properties can be proven from a direct computation using the 
definitions of tab, ^^fc' Qabc ^n.'i tabcd given above but this results in involved 
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calculations even when done by computer. A simpler procedure is to start with (j6.ip 
and derive the relations 

W = Babcdn''n''n^n'', Va = -BptcdhPyn^n'' , (6.2) 

tab = Bpgcdn^n'^hP^h''^, i*^ = BrpsqU^n'^h'' ah\, Qabc = -Bpqrsn^h'^ah''b^\i (6-3) 

tabcd = BpgrshPah\h-,h'a. (6.4) 
From these relations and the properties of the Bel tensor it is straightforward to prove 
the proposition. D 

Remark 6.1 An important consequence of the algebraic properties presented in this 
last result is that tabcd, Qabc t*af^ contain all the information about the Bel 
tensor. As we did in the case of the Bel- Robinson tensor we can count the number of 
independent components of these tensors and check that they add up to the number 
of independent components of the Bel tensor 

number of independent components of Babcd = 45, 

number of independent components of tabcd = 21, 

number of independent components of Q^bc = 18, 

number of independent components of = 6. 

21 + 18 + 6 = 45. 

Proposition 6.2 

W = Q -4=^ tab — ■4=4' tabcd = ■4=4' Rabcd — 0, {no superenergy <J=^ no gravitation) 



Proof : If Rabcd vanishes then so does Babcd and trivially W = 0, tab — 0, tabcd — 0. 
Assume now that W is zero. In that case point plt| of theorem 14. II entails Rabcd = 
thus proving the desired result. D 
We finish this section by pointing out that whenever the Bel and Bel-Robinson 
tensors are equal then we deduce the relations 

tab — tab — ^afci Qabc — Qabc: tabcd — tabcd-; 

from which we conclude that W = W, Va = Va- The Bel and the Bel- Robinson 
tensors are equal if and only if Rab = (see Corollary 6.1 of [38]). 



7. Dynamical lav^rs of superenergy 

In this section we present the most important result of this paper which is the 
orthogonal splitting of (|4.17p . As explained before this result is analogous to (|3.7p - 
(|3.9p and this analogy will enable us to extract some interesting conclusions as to the 
interpretation of certain parts of the orthogonal splitting of the Bel tensor. 

Before presenting the results we should make some remarks concerning the 
calculations. In order to work out the orthogonal splitting of (|4.17p neither (|4.ip . 
nor its orthogonal splitting is needed. This is similar to electromagnetism, where 
the Maxwell equations are not needed to obtain p.7p - (|3.9p . The orthogonal splitting 
of (|4.17p is calculated by inserting the orthogonal splitting of each of the quantities 
appearing in this equation (the Bel tensor, the Riemann tensor and the matter current) 
and then using the orthogonal splitting of the different terms which appear in the 
resulting expressions. Here we only provide the final expressions referring the reader 
to appendix A for more details about the intermediate steps in the calculations. 
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Theorem 7.1 (Dynamical laws of superenergy) The equation 

^a^%cd = ^d'^Rbeca +Zc^Rbeda — -jgcdZ'^^^ Rbfae 

is equivalent to the following set of expressions 

£nPc + Dat*\ + {v. + e - j^Yae + 2^aeC ' Z^e^caP""" + 

+2t^^^ (Q,,, - KJP,) + {h,JV + + 2tta) = 0, (7.1) 

{h^.i-X'^aKe + ^6e + Ybe) ' KeXbc) , (7-2) 

cd + '^''^^afcd + (^"'^'^afcd + -^{tcd + tlj)9 + 2A''{ha(d'Pc) + Q acd) - 
-j^"" (l^afhhcd + 2/l/(d£e)a/.) Z% + {h^dY^e ' /ida^e - hM = 0, (7.3) 

£ntU + -DaQ^fed + '^''^Kfbd + ^"^KfM + I (hd + nd + tbda) + 

{^^a(d'Pb) + Qabd + 0(6d)a) ^ .f'^'^£ha{dhb)fZ°'c " L°-£ae(dZ%-) + J"''' ha(dYb)e + 
+ {ha{dXb)e + he(^d{Xh)a ~ ^b)a-'^^p) " hbdi^ae — haeX^p) — haeXbd) J =0, (7.4) 
£nQbcd + ^at°^bcd + '^"'^^afbcd + '^"'^ ^afbcd + (^Q{bcd) + ^"(^abcd + habtcd + ^haidt*c)b) + 
+rfZ\ Qs^e/Z^ccJ + 2/l,(<i£,)„/) - T'sbefihcdZi - 2h,^dZ{^) + 

+L''{h,dYab - 2/l„(dFe)6) + f'^Haefbcd = 0, (7.5) 

where 

Sa/cd = -2/l^(^tc)/ + '2ha{dtc)f + icdaf, i^afcd = -2h^i^^c)f " '^K(dt*c)f, (7.6) 
^*afbd = tabdf + ha(d(tb)f - K^f), ^Ifbd = -ha{d(ib)f + Sife)/), 

Ha/bcd = -'^Qfcd^ab " (bf)(chd)a-< ^afbcd = 4:Q[fb](chd)aJ 
Haefbcd = 2habhe(dXc)f + 2hejha{dXc)b + 2hadhc[bX f]e + 4ft.c[d/ia] [/-'^6]e + 
+X\{habhcdhef - 2habhc{fhe)d " '^hefha{dhc)b) 

Proof : See appendix A. D 

Theorem 7.2 (Matter current conservation) The equation Va^bc" = Q is 
equivalent to the expressions 

£nLp = A«(Jpg - Jpg) + -(j/g - 2Lp) + {Lbhpa - japb)(^"^ + [japb + Lbhpa)uj"-^ - 

-D,Jp\ (7.7) 
£nJbp = -Dgifc/ + - (2J[bp] - Jbp^ e + Ai {2hg^Lb] - jbpq) + 

+2a°-'' (h^^pJb]c + ha\pJb]c^ + 2ui°-'' {^h^^Jb]c + ha[pJb]^ ■ (7.8) 
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Proof : Again see appendix A. D 

Remark 7.1 Equations ([711) - (|73|) and (fr71) - (fr5|) can be regarded as the 
gravitational counterpart of (|3.7[) - (|3.9[) . They form an inhomogeneous evolution system 
for the variables Va, tab, ^abi Qabc^ and Jab- The inhomogeneous part (source) of 
each equation consists of those terms which contain neither kinematical quantities nor 
spatial covariant derivatives. These terms play the same role as —EaJ"' in (|3.7p (power 
lost by the charge flux) and SacbB^J'^ + EaP (change of momentum due to charges) 
in (|3.8p . We also find that no expressions for £ntabcd, £nJab, £njabc are supplied by 
the orthogonal splitting of (|4.17p and in fact only by using the full content of (|4.ip 
can such expressions be found. 

Remark 7.2 The evolution equations of theorems 17.11 and 17.21 are written in such 
a way that the coupling of the kinematical quantities to the different parts of the 
orthogonal decomposition of the Bel tensor and the matter current is manifest. Note 
also that in these equations we can find terms which do not contain kinematical 
quantities. As the kinematical quantities can be always set to zero at a a given point 
by choosing a suitable vector field n° we deduce that any term containing explicitly a 
kinematic quantity is observer dependent and it will play a similar role as the inertial 
terms in equations (|3.7p - (|3.9p found for electromagnetism. 

Taking the trace of (|7.3p we find 

i:„W+i?a^Va'^'=(iae+2t:j + y(t*'',+2T^)+-£„/b/t^'^^+r''V,e+4A''Pa = 0.(7.9) 

Equations similar to this one have been used in different places of the literature 
principally with the aim of controlling the evolution of the scalar W [TJ [23 • 

7.1. Dynamical laws of superenergy in vacuum 

Theorem 17.11 assumes a far more simpler form in vacuum because the covariant 
divergence of the Bel tensor takes the simpler form Vay^^d = 0. The specific result 
in this case is given in the next theorem. 

Theorem 7.3 The equation 

^a.T%^j^ — 0, 
is equivalent to the following set of expressions 

£ntcd = -2A''{ha(dVc) + Qcda) + 4t^"'/l„(A)6 ^ \icdO " tcdae^l"' " DaQ\a^ (7.10) 
£nQbcd = "^"'{tbcda + '^h.a{d^bc)) + (i^'"'ha(dQbc)e - (^Qbcd - Darted- (^-H) 

Proof : This can be regarded as a particular case of theorem 17.11 with 2abc = and 

Babcd = Tabcd- This eutails tab — tafj, Qabc = Qabc, tabcd = tabcd, La — 0, J ab = 0, 

Jab = 0, jabc — which used in (|7.3p and (|7.5p leads to (|7.10p and (|7.1ip . Equation 
(|7.2p becomes an identity and (|7.ip is now obtained by taking the trace of (|7.3p . D 
In the particular case studied in theorem (|7.3p we find that (|7.9p and (|7.ip acquire 
simpler expressions which are 

£,,W = -AA'^Va - 2We - 3i''V,e - DaV% (7.12) 

59 

£nVd = -A^iStda + hdaW) - —Vd - ^Qdaea"" + 2V"uJda - Daf^- (7.13) 
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The linearized form of (|7.12p was known to Bel [5] and in fact he took this equation as 
the starting point for a definition of a state of intrinsic radiation for the gravitational 
field in vacuum (see subsection 18.11 for further details). The general form of (|7.12p 
was derived in [32] ■ It is interesting to note the formal analogy of (|7.12p - (|7.13p with 
(|3.7p - (|3.8|) where W and take the role of the electromagnetic energy density and 
the Poynting vector respectively. Although (|7.13p has, as far as we know, never been 
obtained in its complete form, the knowledge of (|7.12p . even in its linearized form, 
shown in equation (jS.ip . has been enough to construct the analogy just mentioned 
and a lot of work has been devoted to studying the behaviour of gravitational systems 
by studying the super-energy density and the super -Poynting vector in the system 
-see for example [25l [131 SSI US]- The results obtained are very suggestive but we 
must note that (j7.12p - (|7.13p are not equivalent to ()7.10p - ()7.1ip which in fact contain 
more information. Therefore, if we are to study gravitational radiation by means of 
techniques involving the study of the evolution of the different spatial parts of Bel- 
Robinson tensor then we should start with the general equations (I7.10p - (|7.1ip . This 
matter is addressed in section [8| 

8. Application: superenergy radiative states of the gravitational field 

In electromagnetism, we speak of electromagnetic radiation to mean that 
electromagnetic energy is traveling from one part of a system to another which in 
turn implies the existence of a flux of energy-momentum. By the Poynting theorem 
this flux is represented by the Poynting vector and thus whenever the Poynting vector 
is not zero at a point we say that electromagnetic radiation is going through that point. 
This statement is observer dependent because in order to define the Poynting vector 
an observer is needed (see equation (j3.6l) V Therefore we may find for example, that 
the Poynting vector is zero for one observer whereas another observer measures a non- 
vanishing Poynting vector. However, there are configurations in which any observer 
will measure a non-vanishing Poynting vector and in these cases it is said that the 
electromagnetic field is in a radiation state at the point. From an algebraic point of 
view this can only happen if the electromagnetic field Fab is singular or null which 
means that it can be written as the exterior product of a null and a spatial vector, 
(see e.g. |31]). 

If we try to follow the same procedure to define gravitational radiation in General 
Relativity we are immediately confronted with the fact that, due to the equivalence 
principle, we can always find an observer who measures no "gravitational energy 
density" at a point, for any quantity with dimensions of energy constructed from 
the metric tensor gab (typically this involves expressions which are quadratic in the 
first derivatives of the metric tensor). This means that in General Relativity we cannot 
pursue the same procedure used to define radiating fields as in electromagnetism if we 
insist upon using quantities with dimensions of energy for this purpose. Of 
course, this does not imply that "gravitational energy" is meaningless and in fact we 
can construct quasilocal and global quantities with dimensions of energy which tell us 
when a gravitational system is radiating. This has been performed for the important 
case of isolated systems where the quantity is the Bondi mass [H |37l [S^ . 

If instead of energy, we use superenergy as a replacement, then the afore- 
mentioned problem disappears and one can use the same ideas as in electromagnetism 
to define radiating gravitational fields or radiating spacetimes in a local way. This 
approach was pioneered by Bel many years ago in [5] and, indeed, the results presented 
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in this section can be regarded as a continuation of Bel's work. We must bear in mind 
all the time that radiating gravitational fields defined in terms of superenergy are in 
principle different from radiating fields defined by a quasilocal energy prescription. To 
find the precise relation between both concepts is an interesting open question which is 
a particular case of a more general problem, namely, the possible relationship between 
superenergy and energy. This is a long standing question which has been already 
largely researched O El [24l [28l [29] (a fuller list of references about this subject can 
be found in [38]). 

8.1. Superenergy radiative states for vacuum spacetimes 

Let us start by reviewing Bel's work about the definition of a radiative spacetime. 
The starting point of Bel's study was the linearized form of (|7.12p . To obtain this 
form, we define a coordinate chart (t, x*), i = 1,2,3 in such a way that d/dt is 
the unit timelike vector n"" and {d/dx^} are spacelike Vi. Next we approximate 
the spatial covariant derivative by a covariant derivative compatible with the frame 
{d/dx^,d/dx^,d/dx^}, and ignore terms containing kinematical quantities. Under 
this approximation, equations (j7.12p - (j7.13p become 



These equations can always be obtained at a given point p of the spacetime if we choose 
an observer n° such that all its kinematical quantities vanish at p (such an observer 
always exists according to the equivalence principle). The first equation of (|8.ip has 
the form of a typical conservation law. The vector is, according to this equation, the 
flux of W (superenergy flux) and whenever is zero we see that W does not change 
for the observer d/dt. According to proposition 15 . 21 the superenergy density W is zero 
if and only if Cabcd vanishes as well and besides W is always nonnegative. Therefore, 
it is possible to take W a,s a replacement for the missing concept of "energy density" 
of the gravitation and we may consider that the existence of a flux of superenergy for 
any observer is an indication of the intrinsic presence of gravitational radiation. These 
ideas led Bel to the following definition ^5j. 

Definition 8.1 (State of intrinsic gravitational radiation, Bel 1962.) We say 

that there is a state of intrinsic gravitational radiation at a point p €z V of a vacuum 
spacetime if Va = Va (ji) does not vanish at p for any . 

A well-known consequence of definition lS.ll is that Petrov types N, II and III are always 
radiative. To show this it is enough to recall that the condition V"' = entails (|5.5p 
which can only be true for either type I or type D. Note that definition [ST] does not 
say anything about the radiative character of Petrov types I and D and in fact a more 
general definition would be needed to decide the issue. To obtain a generalization of 
definition lS.ll is our next task. 

To generalize definition 18.11 we need to use the full information coming from the 
orthogonal splitting of Var*^^^ — and not just (|7.12p which only contains part of 
this information. Theorem l7.3l contains all that is needed in our endeavour. If we wish 
to use the variation of superenergy as a tool to define radiative states then we need to 
find the evolution of a spatial tensor whose vanishing is equivalent to the absence of 
a gravitational field (in vacuum this is just the condition Cabcd = 0). Bel's definition 
is based on the scalar W but proposition (|5.2p tells us that the tensor tab plays a 




3 



3 



(8.1) 



i=l 
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similar role (and besides W is not independent of tat)- The propagation of tab is given 
by (|7.10p and we see that the only term in this equation not affected by kinematical 
quantities (and hence intrinsic) is DaQ\^. 

Definition 8.2 (Intrinsic superenergy radiative state in vacuum) In a vac- 
uum spacetime there exists an intrinsic superenergy radiative state at a point p Cz V if 
Qabci'Ti') does not vanish at p for any unit timelike normal n". 

Remark 8.1 We use the name superenergy radiative state instead of Bel's original 
name of radiative state in order to stress the fact that our definition is based on 
gravitational superenergy. 

Note that there are more tensors which have the relevant properties of tab explained 
above and therefore we could use their propagation as the starting point for a definition 
of superenergy radiative state. The consequence of this is that definition 18.21 admits 
alternative but equivalent formulations. To see an example, consider the spatial tensor 

Wab = EacE\ + BacB\. (8.2) 

Clearly, W^a — ^ ^-i^d Wab = <S=^ Cabcd = 0. Moreover, for any spatial vector x"" , 
Wabx'^x^ is non-negative by inspection. We find that in terms of Wab equation (|7.10p 
takes the equivalent form 

£nWcd = A^i-KdVa + ^hadVc + ^acVd + 2Scda) ' 4.6 ab(dW Uj'' - ^0Wcd + 

+ (h^ + h^ahdbW + 3h,dWab^ + DaS^f, (8.3) 

where 

Scda = '2Bh(dEc)e£a^'' ■ (8-4) 

In view of (|8.3p we deduce that definition 18.21 can be formulated by replacing Qabc 
with Sabc- In fact from (|8.4p and (|5.2p we deduce 

Scda — "^{Qacd ^cdQ fta)' Qacd — '^{Scda hcdS (,(j), 

from which we conclude that both Sabc and Qabc contain the same information and 
thus they should be deemed equivalent. We may expect that any reasonable definition 
of a superenergy radiative state should be formulated in terms of a spatial tensor which 
is equivalent to Qabc- Any such tensor can be regarded as the gravitational equivalent 
of electromagnetism's Poynting vector. The tensor Sabc seems to be the simplest choice 
and one may adopt it as the basic geometric object measuring "superenergy flux" . 

Another interesting aspect of (I7.10p - (|7.1ip or (|8.3p . already pointed out in remark 
17. 2[ is the fact that they are written in such a way that the couplings of the kinematical 
quantities to the different spatial parts of the decomposition of the Bel-Robinson 
tensor are apparent. In our present context these couplings could be interpreted as 
the effect on the superenergy radiation due to the acceleration, the expansion, the 
shear and the rotation. At this point it is instructive to compare equation (|7.10p 
(or its equivalent (|8.3p ) with its electromagnetic counterpart which is (|3.7p . In the 
electromagnetic case we realize that the vorticity has no effect whatsoever on the 
radiation of electromagnetic energy whereas it certainly influences the radiation of 
superenergy because (or equivalently ujab) appears explicitly in (I7.10p . 
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8. 2. Superenergy radiative states for general spacetimes 

Using the ideas explained in the previous section we can formulate a definition of an 
intrinsic superenergy radiative state that is similar to definition 18.21 but valid for a 
general spacetime. In this case we need to study the evolution of a spatial quantity 
which is zero if and only if the Riemann tensor vanishes. As stated in proposition 
16.21 the tensor tab has the required properties and hence the terms appearing in the 
evolution equation oitab should enable us to define the concept of an intrinsic radiative 
state. The evolution equation sought is (|7.3p and hence the inspection of this equation 
leads us to the following 

Definition 8.3 (Intrinsic superenergy radiative state in a general spacetime) 

There exists an intrinsic superenergy radiative state at a point p if for any unit 
timelike vector n'^ it is the case that Qatd''^) does not vanish at p. 

Similar considerations as in the case of definition 18.21 apply here. 
9. Conclusions and open issues 

In this work we have obtained the full orthogonal splitting of the Bel tensor and its 
covariant divergence and we have particularized it to the important case of vacuum 
spacetimes where the Bel tensor becomes the Bel- Robinson tensor. This gives rise to 
the dynamical laws of superenergy. The concept of a superenergy radiative state has 
been introduced. The work just presented opens new research lines which we believe 
are worth exploring. Perhaps one of the most interesting issues is a global formulation 
of the dynamical laws of superenergy complementing the local formulation of theorem 
17.11 Such a global formulation would enable us to apply our techniques to realistic 
astrophysical settings such as oscillating stars, rotating bodies or radiating binary 
systems. 

In this paper we have restricted ourselves to the superenergy defined from the 
Riemann and Weyl tensor but one can define tensors representing superenergy from 
a general field resulting in the superenergy tensor oi that field [38] . In this framework 
it is possible to calculate the covariant divergence of a superenergy tensor and obtain 
an expression similar to the first equation in (|4.17p with the Bel tensor replaced by 
a suitable superenergy tensor. The orthogonal splitting of such an equation would 
yield the dynamical laws of the superenergy associated with that particular field. 
An interesting example concerns the electromagnetic field. In this case a possible 
superenergy tensor is the Chevreton tensor which was first introduced in |15j and 
recently stimulating results about its symmetries and the covariant divergence of its 
trace have been obtained [8]. The Chevreton tensor, like the Bel- Robinson tensor, is 
a rank-four tensor and its covariant divergence couples the Weyl tensor with terms 
which contain covariant derivatives of the Faraday tensor [18] . This suggests a possible 
exchange between the gravitational and the electromagnetic superenergies [351 1301 dS] ■ 
The orthogonal splitting of the covariant divergence of the Chevreton tensor might 
shed light on the nature of this exchange. 

Another important issue is the possible relationship between superenergy and any 
of the available quasilocal concepts of gravitational energy which have been developed 
over the years. This is a topic which has been extensively researched in the past 
and no clear conclusion has been reached. In this work no attempt has been made 
in this direction and our point of view has been to regard superenergy as a physical 
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quantity on its own right. We believe that this idea can be put to work by means of 
the results of theorem 17.11 which would demand a formulation of the dynamical laws 
of superenergy tailored for each physical system under study. However, a relation 
between superenergy and gravitational energy cannot be ruled out and the orthogonal 
splitting of the Bel tensor might bring a new point of view to this old problem. 
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Appendix A. Technical details about the computations 

In this appendix we supply details about the calculations required in this work. In 
order to do so we need to explain some implementation aspects of the system xAct. We 
will limit ourselves to only those issues which are needed in our calculations referring 
the interested reader to p33 | for a full documentation and tutorials about xAct. 

Orthogonal splittings play an essential part in our work and the implementation 
of xAct in regard to this matter is completely adapted to our requirements. The basic 
elements of the orthogonal splitting are defined through the command 

In []:= Def Metric [1, h[-a, -b] , cd, {"l","D"}, InducedFrom->-[g, n}, PrintAs->"h"] 

Here h[-a,-b] represents the spatial metric hab which is constructed from the 
spacetime metric gab (represented in the system by g[-a,-b]) and the unit normal 
vector (represented by n[a]). The operator cd[-a] is the Cattaneo operator Da 
associated with hab- The system is able to handle all the properties of the Cattaneo 
operator explained in subsection l2.2l in a natural fashion. 

The general expression for the orthogonal splitting of any tensor is equation (j2.3p . 
This result is implemented in xAct by means of the command 

In[]:= InducedDecomposition[expr, {h.n}] , 

where expr represents any tensorial expression. The output of InducedDecomposition 
is the result of applying formula (j2.3p to expr. The orthogonal projector operator 
Ph which appears in (j2.3p is also implemented in xAct by means of the command 
Projectorh[expr] where again expr represents an arbitrary tensor. The basic com- 
mands just explained enable us to find efficiently orthogonal splittings similar to equa- 
tion (|2.12p with La replaced by any spatial tensor of higher rank. 

Proof of theorems \7.1\ and \7.S\ 

To prove theorems 17.11 and 17.21 we need to find the orthogonal decomposition of the 
equations shown in (j4.17p . The first step is to replace Babcd, Rabcd and Zabc with 
their orthogonal splittings, eqs. (|6.ip . (14. 6p and (I4.10p respectively. The covariant 
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derivatives of are decomposed according to (|2.4p and VaCbcd is decomposed by 
means of the formula 

After doing these replacements we obtain expressions which contain VaM^, VaT't, 
Vaibc, Vat^^, Vaibcde, aLb, ^ aJbc, ^aJbc, ^ ajbcd- Thcsc are further decomposed 
by following the procedure explained in subsection 12.21 For example, the orthogonal 
decomposition of V aLb is just equation (|2.12p which also holds if we replace Lb with 
Vb- Other orthogonal decompositions needed are 

V Jafa = -'2A'^ncn(bta)d + '^ni^a (\'tb)cd + tb){crcd + ^cd) ) + Dctab + 



+nc (^^tab + 2t(f{aa)d + ^a)d) ~ £nta>}j , (A.l) 

which is also valid if we replace tab by any symmetric spatial tensor and 

Vaifac = 2A'^na{Jd[bn'c]) + 2n[c (^^Jb]ad + Jb\{^ad + t^ad)^ + DaJbc + 

+«a i^'^bcO + 2J^^{ac]d + t^cjrf) - £nJb(^ , (A.2) 

which is true if we replace Jbc with any antisymmetric tensor. The expressions for 
the orthogonal splitting of the remaining covariant derivatives are very long and we 
omit them. Inserting the orthogonal splittings in (|4.17p and rearranging the equations 
obtained as polynomials in n° we obtain the expressions 

An'^n'^n'^ + Bln^n"^ + s'^^n'^^n^ + Cfn^ + cl^^n"^^ + E'""^ = 0, (A.3) 

Cf n^l + I^" = 0, (A.4) 

where all the tensor coefficients of these polynomials are spatial. This implies that 
the coefficients of the polynomials must vanish and these conditions lead us to 

^ = 0, Bl = = 0, C['"^'^ = CJ'^ ^ 0, C^^" = 0, E^""^ ^Q,G\^ 0, /I''"! = = 0.(A.5) 

After some manipulations we find that the condition >1 = is equivalent to (|7.9p . 

Bl^ are equivalent to ^^-^J^, Cf = is equivalent to ([7^, = 
is equivalent to (TTi)) . E'"''^ = is equivalent to (|73)) . = is equivalent to (fTT)) 
and I^'^ = is equivalent to (|7.8p . The condition A = is redundant because it can 
be obtained as the trace of Ci'^ = and therefore we do not need to consider it. We 
have thus recovered all the expressions given in the statements of theorems 17.11 and 
17.21 Note that the polynomials (|A.3p - (jA.4[) are equivalent to each of the equations 
presented in (|4.17p and so is the set of conditions stemming from (jA.5p . D 



Appendix B. Canonical forms for the electric and magnetic parts of Weyl 
tensor in the different Petrov types. 

We present next the canonical forms of the electric and magnetic parts of Weyl tensor 
for the different Petrov types. We follow [5] in our presentation (see also [1^ for 
an equivalent representation of the canonical forms). All the canonical forms are 
written with respect to a certain orthonormal frame O = {e", 62, 63} of spatial vectors 
(canonical frame). 
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Petrov type I 

In this type E^b and Bab take the following form in the canonical frame O 

Eab = diag(^ii, i;22, Eas), Bab = diag(Bii, B22, B33), (B.l) 

with the additional conditions 

Ell + E22 + E33 = 0, Bii + B22 + B33 = 0. (B-2) 

Petrov type D 

This type arises if we set —^Eu = E22 = E33, —^Bu = B22 = B33 in the previous 
case. 

Petrov type II 

The canonical forms for £"06) Bab in the frame O are 

Eab = 



Ell 











-^+^23 


-E'23 





E23 


-^--823 


Bii 











^ - -^23 


B23 





-B23 


-^+E23 



Bab= ^ -^-E23 B23 I . (B.3) 

Petrov type III 

The canonical forms for Eab^ Bab in the frame O are 

^12 

(B.4) 



Petrov type N 

The canonical forms for Eab, Bab in the frame O are 

00 0\ /OOO 
Kb = I E22 -B22 , Bab = B22 E22 I . (B.5) 

—B22 —E22 I \ £^22 ~B22 
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